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Abstract. Phonon dynamics, charge response and the phonon density of states are 
calculated for the high-temperature superconductors HgBa2Cu04 and Bi 2 Sr 2 Cu06 
within a microscopic model for the electronic density response. The results are 
compared with previous calculations for La2Cu04 and Nd2CuC>4. Our main focus 
is on the phononanomalies which are connected with the high-frequency oxygen bond- 
stretching modes (OBSM) found before in our calculations for p-doped La2Cu04 and 
n-doped Nd2Cu04. We investigate the question if the characteristic softening of the 
OBSM and the related strong coupling to the electrons is also present in HgBa2Cu04 
and Bi2Sr2Cu06- In particular the importance of the contribution of the more 
delocalized Cu4s state besides the localized Cu3<i state on the softening is investigated 
and the different anticrossing behaviour due to the presence of several phonon modes 
with the same symmetry as the OBSM is studied. This makes the identification of 
the anomalies quite complicate. Furthermore, the influence of electronic polarization 
processes at ions out of the CuO plane in the ionic layers on the phonon dynamics is 
calculated. In this context the qualitative type of charge response, i.e. the presence 
or not of possible metallic charge fluctuations at the Hg or Bi ion, respectively, linked 
via the apex oxygen to the CuO plane, proves to be very sensitive for certain phonon 
modes. All the calculations are compared to the experimental results available so far. 
The latter, however, are rather incomplete for the Hg and the Bi compound. 
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1. Introduction 

In recent literature there is increasing evidence that the electron-phonon interaction 
(EPI) is significant in the cuprate based high-temperature superconductors (HTSC's) 
and it is argued that phonons might play an important role for the electron dynamics 
in the HTSC's; see, e.g. Refs. [HEIEIHEIEIIZIIS]. Itis known from the 
experiments for some time that the frequencies of the oxygen bond-stretching modes 
(OBSM) are strongly renormalized (softened) upon doping in the p-doped cuprates 
|9] QUI [TTJ [121 EES HH EE] • Quite recently first evidence of an anomalous dispersion of 
the OBSM also has been reported in a Bi cuprate [16J which, among other things, will be 
investigated theoretically in this work. Anomalous softening of the OBSM, not existing 
in the undoped insulating state, also has been observed in the n-doped metallic state 
of Nd2Cu04 [171 [THl HE]. The phonon dispersion and the charge response in Nd2Cu04 
have been calculated in [20] and the result compares well with the measured anomalous 
dispersion. 

All these findings support a generic nature of the OBSM phonon anomalies in the 
HTSC's. The frequency renormalization and a corresponding increase of the linewidths 
point to a strong coupling of these phonons to the charge carriers. Deviations from the 
presumably generic dispersion for the OBSM in form of a sharp softening have recently 
been reported in inelastic neutron scattering (INS) and high resolution inelastic x-ray 
scattering (IXS) experiments, respectively, for La 2 _ x Ba x Cu04 at x — 0.125 [6] and 
x = 0.14 [21]. This has been related to the possibility of stripe formation in this 
compound. 

According to our calculations the origin of the phonon anomalies is a strong nonlocal 
EPI due to specific screening effects in form of a polarization created by metallic charge 
fluctuations (CF's) on the outer shells of the Cu and O xy ions in the CuO plane; see, e.g. 
[221 1231 121] . This type of screening characteristic for a strongly inhomogeneous electronic 
system like the cuprates with a strong component of ionic binding and a corresponding 
density distribution which is mainly concentrated at the ions in the CuO plane leads for 
the p-doped cuprates [221 121] and also for the n-doped Nd 2 Cu0 4 [20] to an anomalous 
behaviour of the dispersion. 

The specific form of screening yields for the OBSM in metallic La2Cu04 and 
Nd2Cu04 a downward dispersion of the anomalous longitudinal branch with high 
frequencies along the (1,0,0) and the (1,1,0) direction, respectively. In case of Nd2Cu04, 
however, the dispersion of the branches with high frequencies comprising the anomalous 
OBSM is more complex than in La2Cu04 because there are in contrast to La2Cu04 
several branches of the same symmetry type nearby in energy leading to complex 
anticrossing phenomena [20]. Such a possibility is also investigated for HgBa 2 Cu0 4 
and Bi 2 Sr 2 Cu06 in this work. 

The anomalous dispersion cannot be understood within a typical lattice dynamical 
model, like the shell model, extended by a homogeneous electron gas screening that is 
diagonal in reciprocal space. This type of dielectric response always yields an increasing 
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dispersion when passing from the center of the Brillouin zone (BZ) into the zone, because 
of an incomplete screening of the phonon- induced changes of the Coulomb potential at 
shorter distances. In contrast, the anomalous softening of the OBSM in particular of 
the oxygen half-breathing mode (Ax/2 anomaly halfway the A ~ (1,0,0) direction) 
and the full oxygen breathing mode (Og at the X point of the BZ) is a result of 
an " overscreening" of the phonon induced changes of the Coulomb potential. The 
overscreening effect has been shown by our calculations to be due to nonlocally excited 
CF's localized at the Cu and O xy sites in the CuO plane accompanied by (dynamic) 
charge ordering in form of localized stripes of alternating sign in the plane [22[ 123] . 
Expressed in terms of the density response matrix in reciprocal space such a strong 
nonlocal EPI effect and "overscreening" of the OBSM is related to the off- diagonal 
elements (local field effect) being apparently very important in the cuprates as revealed 
by the anomalous softening found in our studies. 

In our calculations a sufficiently broad set of orbital degrees of freedom (Cu3d, 
Cu4s, 02p), the full three dimensional long range Coulomb interaction as well as the 
short range local repulsion of the electrons, especially the important on-site repulsion 
mediated by the localized Cu3<i orbital, is considered quantitatively. All the resulting 
couplings arising in the dynamical matrix are microscopically well defined and can be 
calculated. 

From our investigations, see, e.g. Ref. [26J, we extract a strong EPI for the 
phonon anomalies Ax/2 and Og , which can well be treated in adiabatic approximation. 
Moreover, we find a very strong nonadiabatic enhancement of the coupling compared 
to the adiabatic case for phonon-like modes, like Of (axial oxygen breathing mode at 
the Z point of the BZ) , from a small nonadiabatic region around the c-axis. In this 
region poor dynamical screening due to the slow charge dynamics around the c-axis of 
the long range polar Coulomb interaction of the ions leads in metallic LaCu04 to low- 
energy plasmons mixing with certain interlayer phonons of the same symmetry. For a 
discussion of phonon-plasmon mixing in the cuprates, see also [27]. An essential nonlocal 
interaction of holes with c-axis polarized optical phonons in the HTSC's has also been 
reported in [281129] . 

Our calculations show that besides the more delocalized Cu4s state the short range 
part of the Coulomb interaction Ud related to the localized Cu3<i orbital is particular 
important for the Ax/2 and Og mode polarized and propagating in the CuO plane while 
the strong long range nonlocal polar Coulomb interaction plays an important role for Of 
and similar A x phonons polarized and propagating along the c-axis. Thus, the relevance, 
of both, short range electron repulsion Ud in the correlated Cu3<i state as well as long 
range polar coupling of the electrons to the phonons, is demonstrated by our calculations 
of the phonon anomalies in the cuprates. Consequently, any quantitative study of the 
physics of the normal or superconducting state of the cuprates should consider electron 
correlations and strong nonlocal EPI on an equal footing. 

The paper is organized as follows. In section 2 the theory and modeling is shortly 
reviewed. Section 3 presents the calculations for the phonon dynamics, charge response 
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and the phonon density of states of HgBa 2 Cu04 and Bi 2 Sr 2 Cu06 and a comparison 
with earlier calculations of LaCu04 and Nd 2 Cu04. The main focus is on the OBSM 
phonon anomalies in the various materials and the effect of the charge response in the 
ionic layers on the dispersion. A summary of the paper is given in section 4 and the 
conclusions are drawn. 

2. Theory and modeling 

In the following a brief survey of the theory and modeling is presented. A detailed 
description can be found in [23] and in particular in [30] where the calculation of the 
coupling parameters of the theory is presented. 

The local part of the electronic charge response and the EPI is approximated in the 
spirit of the quasi-ion approach [201 EI] by an ab initio rigid ion model (RIM) taking into 
account covalent ion softening in terms of (static) effective ionic charges calculated from 
a tight-binding analysis. The tight-binding analysis supplies these charges as extracted 
from the orbital occupation numbers of the \i (tight-binding) orbital in question: 



C Mn (k) stands for the /^-component of the eigenvector of band n at the wavevector k in 
the first BZ; the summation in equation ([T]) runs over all occupied states and N gives 
the number of the elementary cells in the (periodic) crystal. 

In addition, scaling of the short range part of certain pair potentials between the 
ions is performed to simulate further covalence effects in the calculation in such a way 
that the energy-minimized structure is as close as possible to the experimental one 
[32]. Structure optimization and energy minimization is very important for a reliable 
calculation of the phonon dynamics through the dynamical matrix. Taking just the 
experimental structure data as is done in many cases in the literature may lead to 
uncontrolled errors in the phonon calculations. 

The RIM with the corrections just mentioned then serves as an unbiased reference 
system for the description of the HTSC's and can be considered as a first approximation 
for the insulating state of these compounds because of the strong ionic nature of 
bonding in the cuprates. Starting with such an unprejudiced rigid reference system non- 
rigid electronic polarization processes are introduced in form of more or less localized 
electronic charge-fluctuations (CF's) at the outer shells of the ions. Especially in the 
metallic state of the HTSC's the latter dominate the nonlocal contribution of the 
electronic density response and the EPI and are particularly important in the CuO 
planes. In addition, anisotropic dipole-fluctuations (DF's) are admitted in our approach 
[30], [33] . They prove to be specifically of interest for the ions in the ionic layers mediating 
the dielectric coupling and for the polar modes. Thus, the basic variable of our model 
is the ionic density which is given in the perturbed state by 




(1) 



rik 



P«(r, Qx, Pa) = p° (r) + QxPx F (r) + p a • rp°(r). 



(2) 



A 
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p° a is the density of the unperturbed ion, as used in the RIM, localized at the sublattice 
a of the crystal and moving rigidly with the latter under displacement. The Q\ and p^ F 
describe the amplitudes and the form-factors of the CF's and the last term in equation 
(J2]) represents the dipolar deformation of an ion a with amplitude (dipole moment) p a 
and a radial density distribution p®. r denotes the unit vector in the direction of r. The 
are approximated by a spherical average of the orbital densities of the ionic shells 
calculated in LDA taking self- interaction effects (SIC) into account. The dipole density 
is obtained from a modified Sternheimer method in the framework of LDA-SIC [30J. 
All SIC-calculations are performed for the average spherical shell in the orbital-averaged 
form according to Perdew and Zunger [31]. For the correlation part of the energy per 
electron the parametrization given in [34J has been used. 

The total energy of the crystal is obtained by assuming that the density can be 
approximated by a superposition of overlapping densities p a . The p° a in equation (j2J) 
are also calculated within LDA-SIC taking environment effects, via a Watson sphere 
potential and the calculated static effective charges of the ions into account. The 
Watson sphere method is only used for the oxygen ions and the depth of the Watson 
sphere potential is set as the Madelung potential at the corresponding site. Such an 
approximation holds well in the HTSC's [32j [35]. Finally, applying the pair-potential 
approximation we get for the total energy: 

E{R,Q = Y,K(Q + \ E <MBJ-H£,C). (3) 

a,a (a,«)^(b,/3) 

The energy E depends on both the configuration of the ions {R} and the electronic 
(charge) degrees of freedom (EDF) {(} of the charge density, i.e. {Q\} and {p Q } in 
equation <^I§. E* are the energies of the single ions, a, b denote the elementary cells and 
a, (3 the corresponding sublattices. The second term in equation ([3]) is the interaction 
energy of the system, expressed in terms of anisotropic pair- interactions $> a p. Both E* 
and & a/ 3 in general depend upon ( via p a in equation (j2J). 

The pair potentials in equation (jBJ) can be seperated into long range Coulomb 
contributions and short range terms, for details see e.g. Ref [30J. 

From the adiabatic condition 

^ = (4) 

the electronic degrees of freedom ( can be eliminated, an expression for the atomic force 
constants can be given and from this the dynamical matrix in harmonic approximation 
can be derived as 



#(q) 



RIM 



The first term on the right hand side denotes the contribution from the RIM. M a , Mp 
are the masses of the ions and q is a wave vector from the first BZ. 
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The quantities B(q) and C(q) in equation ^ represent the Fourier transforms of 
the electronic coupling coefficients and are calculated from the energy in equation fl3]), 
i.e. 



,ab d*E(R,() 



k denotes the EDF (CF and DF in the present model, see equation (J2J)) in an elementary 
cell. The B coefficients describe the coupling between the EDF and the displaced 
ions (bare electron-phonon coupling), and the coefficients C determine the interaction 
between the EDF. The phonon frequencies u a (q) and the corresponding eigenvectors 
e Q (q<r) of the modes (qcr) are obtained from the secular equation for the dynamical 
matrix in equation dSJ), i.e. 

5Z^/(q)ef(q) =o; 2 (q)er(q). (8) 

The Eqs. flSJ-flB]) are generally valid and, in particular, are independent of the specific 
model for the decomposition of the perturbed density in equation (j2J) and the pair 
approximation in equation ([3]) for the energy. The lenghty details of the calculation of 
the coupling coefficients B and C cannot be reviewed in this paper. They are presented 
in [30]. In this context we note that the coupling matrix C KK >(q) of the EDF-EDF 
interaction, whose inverse appears in equation (JSJ) for the dynamical matrix, can be 
written in matrix notation as 

c = ir 1 + v. (9) 

IT 1 is the inverse of the proper polarization part of the density response matrix and 
contains the kinetic part to the interaction C while V embodies the Hartree and 
exchange-correlation contribution, because they are related to the second functional 
derivatives with respect to the density p of the kinetic energy and the exchange- 
correlation energy, respectively [23]. C^ 1 needed for the dynamical matrix and the 
EPI is closely related to the linear density response matrix and to the inverse dielectric 
matrix e" 1 , respectively. 

Only very few attempts have been made to calculate the phonon dispersion and 
the EPI of the HTSC's using the linear response method in form of density functional 
perturbation theory (DFPT) within LDA [361 E3 EH [39] . These calculations correspond 
to calculating IT and V in DFT-LDA and for the metallic state only. On the other 
hand, in our microscopic modeling DFT-LDA-SIC calculations are performed for the 
various densities in equation (j2J) in order to obtain the coupling coefficients B and V. 
Including SIC is particularly important for localized orbitals like Cu3d in the HTSC's. 
Our theoretical results for the phonon dispersion J22J [2U 123 [33], which compare well 
with the experiments, demonstrate that the approximative calculation of the coupling 
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coefficients in our approach is sufficient, even for the localized Cu3g? states. Written in 
matrix notation we get for the density response matrix the relation 

C- 1 = U(l + VUy 1 = Ue~\ e = l + VTl. (10) 

The CF-CF submatrix of the matrix II can approximately be calculated for the metallic 
(but not for the undoped and underdoped) state of the HTSC's from a TBA of a single 
particle electronic bandstructure. In this case the electronic polarizability IT in tight- 
binding representation reads: 

2 ^ /n'(k+q)-/„(k) 



U KK ,(q,u = 0) = -— 



N ^ R n ,(k + q)-£ n (k) 

n,n ,k 



X 



x [C: n (k)C Kn ,(k + q)] [C£ B (k)CW(k + q)]* . (11) 

/, E and C in equation (iTTj) are the occupation numbers, the single-particle energies and 
the expansion coefficientes of the Bloch-functions in terms of tight-binding functions. 

The self-consistent change of an EDF at an ion induced by a phonon mode (qc) 
with frequency w CT (q) and eigenvector e°(qcr) can be derived in the form 



J>™(qK( q( x) 



= 5Uqo-)e^, (12) 
with the displacement of the ions 

<™ = fcdbr e " (q ' 7)e " ,R * 

= u a ( v y> K \ (13) 

The self-consistent response per unit displacement of the EDF in equation ( I12jl is 
calculated in linear response theory as: 

X(q) = n(q)e- 1 (q)B(q) = C'- 1 (q)B(q). (14) 

The generalization for the quantity II in Eqs. (Q and (|T0|) needed for the kinetic part of 
the charge response in the nonadiabatic regime, where dynamical screening effects must 
be considered, can be achieved by adding {hui+irj) to the differences of the single-particle 
energies in the denominator of the expression for II in equation ( TTTT) . Other possible 
nonadiabatic contributions to C related to dynamical exchange-correlation effects and 
the phonons themselves are beyond the scope of the present approach. Using equation 
( JTOl) for the dielectric matrix, e, and the frequency-dependent version of the irreducible 
polarization part, II, according to equation (ITTj) . the free-plasmon dispersion is obtained 
from the condition, 

det[e we ,(q,w)] = 0. (15) 

The coupled-mode frequencies of the phonons and the plasmons must be determined self- 
consistently from the secular equation (JSJ) for the dynamical matrix which now contains 
the frequency u implicitly via II in the response function C^ 1 . Such a nonadiabatic 
approach is necessary for a description of the interlayer phonons and the charge-response 
within a small region around the c-axis as performed in [221 EE] . 
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3. Results and discussion 

3.1. Ionic reference system and structural data 

For a definitive discussion of the renormalization of the phonon modes introduced by 
the nonlocal EPI effects related to the electronic polarization processes of CF and DF 
type, a quantitative reference model for the calculation of the phonon dispersion based 
essentially on the important component of ionic binding is necessary. A suitable model 
sketched in section 2 which represents approximately the local EPI effects is provided 
by the ah initio rigid-ion model, extended by covalent ion softening and scaling of the 
short range part of certain pair potentials. 

A modified ionic description is quite general a good starting point for the physics in 
the cuprates in particular for their insulating state. In this context we also refer to the 
work in |40j HTJ H2] . In our calculations we found a suitable set of static effective charges 
given in table [1] which together with covalent scaling leads to a good overall agreement 
with the experimental structural data, shown in table [2j The charges of the Cu and the 
O xy ions in the CuO plane are taken identical for La 2 Cu0 4 , Nd 2 Cu0 4 and Bi 2 Sr 2 Cu0 6 . 
They result from a tight-binding analysis of the ab initio band structure in La 2 Cu04 
|32j. This means that the CuO plane is treated as generic for these materials as far as 
the effective charges are concerned. On the other hand, for HgBa 2 Cu04 these charges 
have to be chosen more ionic, otherwise the phonon modes are too low in frequency as 
compared with the experiment. In all the materials covalent scaling of the short range 
pair potential between the Cu and the O xy ions in the plane is very essential and in case 
of HgBa 2 Cu04 and Bi 2 Sr 2 Cu0 6 in addition between the apex oxygen O z and Hg and 
O z and Bi, respectively. 

This is consistent with a significant shorter distance of the 2 -Hg and O z -Bi link 
as compared to O z -L& or O z -Nd in La 2 Cu04 and Nd 2 Cu0 4 , respectively (O z -Hg: 1.97 
A; O z -Bi: 1.94 A; O z -L&: 2.38 A; 2 -Nd: 2.32 A). Thus, the overlap of the orbital 
densities in the O^-Hg and O z -Bi bond is sufficiently large and consequently metallic 
CF's between the partially occupied O z 2p and Hg6s orbitals and the O z 2p and Bi6p 
orbitals, respectively, can be expected. The important effect of such CF's on certain 
phonon modes will be discussed in section 3.2. 



La 2 


Cu0 4 


Nd 2 Cu0 4 


HgBa 2 Cu0 4 


Bi 2 S 


r 2 Cu0 6 


La: 


2.28+ 


Nd: 2.35+ 


Hg: 


1.44+ 


Bi: 


2.28+ 








Ba: 


1.69+ 


Sr: 


1.85+ 


Cu: 


1.22+ 


Cu: 1.22+ 


Cu: 


1.70+ 


Cu: 


1.22+ 


Oxy 


1.42- 


O xy : 1.42- 


Oxy 


1.56- 


Oxy 


: 1.42- 


O z : 


1.47- 


2 : 1.54- 


O z : 


1.70- 


2 : 


1.47- 



03: 1.85- 



Table 1. Static effective charges of the ions. 
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La 2 Cu0 4 p] Nd 2 Cu0 4 [S] HgBa 2 Cu0 4 (45j Bi 2 Sr 2 Cu0 6 [S] 



a = 


3.763 


a = 


A 1 f f 

4.155 


a = 


3.8995 


a = 


3.741 




/o 01 a\ 

(3.810) 




(3.945) 




(3.875) 




(3.798) 


c = 


13.197 


c = 


12.102 


c = 


9.521 


c = 


23.490 




(13.240) 




(12.171) 




(9.513) 




(24.662) 


c/a = 


3.5 


c/a = 


2.9 


c/a = 


O A 

2.4 


c/a = 


6.3 




(3.5) 




(2.9) 




(2.5) 




la n\ 

(6.5) 


t(O z ) = 


0.185 


t(O z ) = 


0.250 


t(O z ) = 


0.313 


t(O z ) = 


0.109 




(0.182) 




(0.250) 




on a \ 

(0.294) 




(0.105) 


r(La) = 


0.134 


/TV T 1 \ 

r(Nd) = 


0.153 


r(Ba) = 


0.183 


r(03) = 


0.199 




(0.138) 




(0.149) 




(0.198) 




/n 1 o/ i \ 

(0.186) 














r(Bi) = 


0.189 
















( C\ A O A\ 

(0.184) 














r(Sr) = 


0.068 
















(0.071) 


C\l-0 X y = 


1.881 


CU-0 X y = 


2.078 


Cu-O xy = 


1.950 


Cu-O x?/ = 


1.870 




(1.895) 




(1.972) 




(1.940) 




(1.896) 


Cu-O z = 


2.441 


Cu-O z = 


3.670 


Cu-O z = 


2.980 


Cu-0 2 = 


2.556 




(2.407) 




(3.622) 




(2.793) 




(2.585) 


O z -La= 


2.382 


2 -Nd= 


2.389 


O z -Hg= 


1.800 


2 -Bi= 


1.893 




(2.381) 




(2.323) 




(1.966) 




(1.944) 



Table 2. Structural data for the investigated compounds. For each lattice parameter, 
the first row shows the calculated data and the second row the experimental data 
in brackets. The lattice constants a and c and the distances between the ions are 
given in units of A, the relative translation t is in units of c. In Nd2CuC>4, t(O z ) is 
fixed by symmetry. The experimental references are listed following the name of the 
compounds. 

3.2. Phonon dispersion, charge response and phonon density of states 

The calculation of the phonon dispersion in this section is based on the dynamical matrix 
given in equation (JSJ) within adiabatic approximation which has been shown to be a good 
approximation for all phonon modes outside a small region around the c-axis, where, on 
the other hand, a nonadiabatic approach is necessary [22J, [26] . As discussed, e.g. in Ref. 
[26J , the nonadiabatic region around the c-axis is so small that it cannot be resolved by 
INS experiments or other probes so far. Thus, only an average of the dispersion can be 
measured along the A ~ (0, 0, 1) direction for modes of Ai symmetry because of their 
coupling to a low lying c-axis plasmon as discussed in [26J. Consequently, all results 
for the Ax-phonons propagating and polarized along the c-axis calculated in adiabatic 
approximation or reported in the experiments must be interpreted as an average over 
the small region with a nonadiabatic charge repsonse. 

Anisotropic DF's being particular important along the c-axis and for ions in the 
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(a)HBM (b)12BM 

Figure 1. Electronic band structure of the (a) 11BM and (b) 12BM. 



ionic layers [201 [33j H7] have been taken into account in all computations of the phonon 
modes in this paper. 

The subsequent calculations are representative for the well doped metallic state of 
the materials where a fully developed Fermi surface exists. Investigations of the charge 
response and phonon dynamics of the underdoped pseudogap state and the insulating 
state of the cuprates have been performed in [201 I2U E3] . These calculations are based 
on a microscopic modeling of the electronic polarizability in terms of rigorous orbital 
specific compressibility-incompressibility sum rules for Il KK /(q) in the long wave-length 
limit in order to discriminate between the charge response of the different electronic 
sates of the cuprates. 

In the following we discuss the effect of the renormalization of the phonon 
dispersion by metallic CF's which are of particular importance for the anomalous 
high-frequency oxygen bond-stretching modes OBSM. In our earlier calculations of the 
phonon anomalies and the phonon dispersion of La 2 Cu0 4 and Nd 2 Cu0 4 [20 1 124 1 [33 | |4"8] 
we have applied successfully a parameterized description (Il-model) for the electronic 
polarizabality matrix I1 KK / in equation (TTTT) . This model highlights besides the effect 
of the localized Cu3g? state the importance of the delocalized Cu4s-component in the 
charge response and for the softening of the anomalous OBSM phonon modes. The 
importance of the Cu4s orbital for a realistic description of the electronic structure of 
the HTSC's has also been pointed out in Refs. [HH [50] . 

In order to demonstrate selectively the renormalization effect introduced by the 
Cu4s orbital for the OBSM not only within the Il-model approach we have calculated 
Ii KK i approximately using for the electronic band structure of the CuO plane an eleven- 
band-model (11BM: based on Cu3<i and 02p states [51]) and a twelve-band-model 
(12BM) which takes the Cu4s state in addition into account. 

In figured] the result for the electronic band structure of the 11BM and the 12BM 
is displayed and in figure [5] the corresponding results for the OBSM phonon anomalies 
are shown for the case of L^CuO^ The important contribution of the Cu4s degree 
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Figure 2. Calculated OBSM phonon anomalies in La 2 Cu0 4 : RIM (— • • — ), RIM 

with anisotropic DF's ( ), 11BM ( ), 12BM ( ) and Il-model ( ). 

The squares (■) represent experimental results [5]. 

of freedom in the charge response for the anomalies is evident by inspection of the 
dotted and broken curve in figure [2j The full curve relies on a diagonal Il-model [48] 
where the matrix elements of the polarizability matrix U(Cu3d) = 2.8 eV _1 , Il(Cu4s) 
= 0.055 eV^ 1 and U(O xy 2p) = 0.2 eV _1 have been extracted from the full polarizability 
matrix, as obtained from a tight-binding analysis of the electronic band structure for 
La 2 Cu04 (31BM [52]) within the framework of the local density approximation (LDA) 
of density functional theory (DFT). The value for Il(Cu4s), however, has been slightly 
increased compared with the 31BM from 0.05 eV" 1 to 0.055 eV _1 in order to optimize 
the calculated result for the A x branch in comparison to the experiment. 

Using the 11BM or the 12BM, respectively, as generic models for the electronic 
structure of the CuO plane also in the calculation of the phonon dispersion of Nd 2 Cu0 4 , 
HgBa 2 Cu04 and Bi 2 Sr 2 Cu06 we find that the softening of the OBSM due to the Cu4s 
CF's is of the same size as for La 2 Cu0 4 displayed in figure [2j 

The detailed results for the renormalization of the OBSM are shown in figure [3] 
where the displacement patterns of the half-breathing mode (Ai/2 anomaly) and the full 
breathing mode Og are given. Moreover, we have listed the corresponding frequencies 
in THz as calculated in the RIM, the 11BM and the 12BM. It should be noted, that 
due to the simple tetragonal elementary cell, in HgBa 2 Cu04 the full breathing mode is 
located at M and the half-breathing mode at X, but to be consistent with the other 
three compounds we keep the denotation A t /2 and Og also for HgBa 2 Cu0 4 . 

A measure for the electron-phonon coupling strength is given by the mode- 
and orbital selective charge fluctuations 5( K (qa) from equation ( fl2l) and the related 
renormalization of the corresponding phonon modes. We find a very strong 
renormalization of the mode frequencies by the electronic polarizations processes which 
nearly exclusively results from the nonlocally excited CF's 5( K listed in table [31 The 
softening for Og as compared to the RIM is in the range between 4 THz and 5 THz and 
for Ai/2 between 2 THz and 4 THz. The largest softening for Ai/2 occurs in La 2 Cu04 
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22.065 16.507 22.977 18.410 

18.988 14.970 20.178 16.647 

17.733 14.371 18.932 16.233 



Figure 3. Displacement patterns of the OBSM for (a) La2Cu04, (b) NC12CUO4, (c) 
HgBa2Cu04 and (d) Bi 2 Sr 2 Cu06. The left pattern shows Og and the right pattern 
Ai/2, respectively. For each displacement pattern the phonon frequency is given for 
the RIM (first row), the 11BM (second row) and the 12BM (third row) in units of 
THz. 

partly due to the fact that there is virtually no anticrossing with other modes of Ai 
symmetry. 

The significant renormalization of the OBSM in all the compounds investigated 
so far signals a strong nonlocal electron-phonon coupling in the cuprates. This strong 
coupling, in particular the effect of the Cu4s orbital, can also be read off from the 
calculated magnitude of the phonon induced charge response given in table |3] in terms 
of the CF's 5( K (qo~) according to equation ffl2"]) . 

Our calculations for HgBa 2 Cu0 4 and Bi 2 Sr 2 Cu0 6 in the following, are based on the 
Il-model approach because unlike to La2Cu04 no reliable tight-binding parametrization 
of a first principles band structure is available for these compounds. Even if a tight- 
binding parametrization of a typical DFT-LDA based band structure would exist it 
would be not sufficient to describe the dynamics of the Ai modes in the HTSC polarized 
and propagating along the c-axis. This is, because typical DFT-LDA calculations are 
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La 2 Cu0 4 

CuM Cu4s O y 2p 


Nd 2 Cu0 4 

Cu3d Cu4s O y 2p 


C) x 

Ax/2 


-23.216 / 
-19 949 -13 255 - 
-16.946 / -6.520 
-15.393 -8.685 -3.283 

HgBa 2 Cu0 4 

Cu3d Cu4s O y 2p 


-18.860 / 
-16 976 -13 129 - 
-14.221 / -5.666 
-13.284 -8.473 -3.152 

Bi 2 Sr 2 Cu0 6 

Cu3d Cu4s O y 2p 


Ax/2 


-24.231 / 
-17.100 -13.576 
-17.528 / -8.865 
-16.276 -9.228 -5.904 


-22.078 / 
-17.081 -11.575 
-13.881 / -8.756 
-11.955 -6.473 -4.002 



Table 3. Magnitudes of the calculated charge CF's <5£ K in units of 10 -3 particles 
for the breathing and half-breathing modes. In the first row the values for the 11BM 
are given, the second row shows the values for the 12BM, respectively. The symbol 
/ means, that CF's are not allowed in the model, — means CF's are not excited for 
symmetry reasons. The corresponding displacement patterns are shown in figure [3l 

much too isotropic and underestimate the anisotropy of the real materials by far. For a 
detailed discussion of this interesting topic in case of La 2 Cu0 4 and its relevance for the 
nonadiabatic charge response and the related phonon-plasmon mixing, see Refs. [26, 53J. 
A short comment should be made concerning the size of the electron-phonon coupling 
strength of the adiabatic OBSM with amplitudes within the plane (Ax/2, O x ) compared 
to modes with out of plane amplitudes as the nonadiabatic apex oxygen bond-stretching 
mode Of. Like in the present paper in our calculations for La 2 Cu0 4 [26] an already 
strong coupling of the adiabatic modes Ax/2 and O x has been found. However, for the 
nonadiabatic Of mode with out of plane amplitudes we obtain in addition a very strong 
increase of the coupling as compared to the adiabatic in plane modes. 

Phonon-plasmon mixing in the nonadiabatic region around the c-axis provides an 
extra channel for the formation of cooper pairs. For a suitable strength of interlayer 
coupling it becomes possible to achieve large contributions to the coupling from both, 
the phonon-like and the plasmon-like modes, as discussed in [221 ES] . Assuming, that 
the phonon-plasmon contribution to pairing is significant this finding could help in the 
search of new high-T c materials. This can be accomplished by changing the out of plane 
solid-state chemistry in order to tune the interlayer coupling. 

Within our H-model approach we can model the charge response for ions out of the 
CuO plane by adjusting the corresponding polarizability matrix elements of these ions 
in order to obtain a reasonable description of the frequency of certain phonon modes 
which are most sensitive to the charge response of the most important out-of-plane ions. 
So, as a first approximation we take the H-model for the CuO plane as introduced for 
La 2 Cu0 4 above also for HgBa 2 Cu0 4 and Bi 2 Sr 2 Cu06- From table H] a comparison of 
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20.009 12.496 20.634 14.083 

17.695 12.363 18.780 12.164 



Figure 4. Displacement patterns of selected optical modes for (a) HgBa2Cu04 and 
(b) r^S^CuOg. The left pattern shows the Raman active mode and the right 
pattern the infrared active Am mode discussed in the text, respectively. For each 
displacement pattern the phonon frequency is given for model II-2dim (first row) and 
model II-3dim (second row) in units of THz. 



HgBa 2 Cu0 4 




M 9 (O r z ) 


Eg 


Eg 


A 2u 


A 2u 


A 2u 


A 2u 




n-2dim 


5.17 


20.01 


2.61 


4.20 


3.44 


4.87 


12.50 


17.55 




Il-3dim 


4.63 


17.70 


2.61 


4.20 


3.31 


4.86 


12.36 


17.32 




Exp. [56] 


4.83 


17.75 


2.25 


5.04 












Exp. [59] 










2.58 


4.50 


10.67 


18.80 




Bi 2 Sr 2 Cu0 6 






Alg 


M 9 (O r z ) 


A 2u 


A 2u 


A 2u 


A 2u 


A 2u 


n-2dim 


4.18 


7.02 


14.06 


20.63 


2.99 


5.09 


9.25 


14.08 


18.67 


H-3dim 


4.16 


5.89 


13.59 


18.78 


2.79 


5.08 


9.17 


12.16 


18.10 


Exp. [57] 


2.07 


3.57 


13.61 


18.71 












Exp. [58] 


3.60 


6.00 


13.76 


18.74 












Exp. [60] 










3.27 


4.95 


9.02 


11.54 


17.57 



Table 4. Frequencies of the calculated and experimental measured Raman- and 
infrared-active modes in units of THz. 



the calculated frequencies of the Raman modes with the experiment shows that the 
largest difference between theory and experiment appears for the Raman mode with 
the highest frequency. Here the apex oxygen ion vibrates along the c-axis against the 
Hg or Bi ion, respectively, which are located just above the apex oxygen, see figure HJ 
As already mentioned in section 3.1 there is a relatively short distance between and 
Hg or Bi, respectively, and consequently metallic CF's can be expected which should 
renormalize the corresponding modes. 

In order to investigate this effect we extend the n-model for the CuO plane (n~ 
2dim) by allowing for metallic charge fluctuations on the O z , Hg and Bi ions (n-3dim). 
The corresponding matrix elements n(Hg6s), U(O z 2p), for HgBa2CuC>4 and n(Bi6p), 
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Figure 5. Calculated phonon dispersion of HgBa2Cu04 for model II-3dim in the main 
symmetry directions A ~ (1,0,0), £ ~ (1,1,0) and A ~ (0,0,1). The symbols o and 
□ represent experimental results from |54) and |12j . respectively. The circle • marks 
two breathing modes at M and the half-breathing mode at X. The arrangement of 
the panels from left to right according to the different irreducible representations is 

as follows:|A 1 |A 2 ( ), A 4 ( )|A 3 |S 1 |S 2 ( ), S 4 ( ) |E 3 |Ai ( ), A 2 

( ) | -A3 1 . Imaginary frequencies are represented as negative numbers. 



U(O z 2p) for Bi 2 Sr 2 Cu06 are determined in such a way that the calculated and measured 
frequencies for are very close, see table ES and figure HI 

Moreover, we extract from the results in table H] that also the other calculated 
Raman modes are improved by allowing CF's at these ions out of the CuO plane. In 
case of Bi 2 Sr 2 Cu06 it is remarkable that also an infrared active A 2u mode is softened 
by nearly 2 THz if the out-of-plane CF's are in addition admitted. As can be expected, 
in this mode the O z and the Bi ion are vibrating against each other, compare with 
figure HI It should be noted that in calculations of the infrared and Raman active 
modes for La 2 Cu0 4 and Nd 2 Cu0 4 with CF's allowed in the Ii-model on O^, La or Nd, 
respectively, these CF only do have a minor influence on the phonon frequencies. 

In figure [5] and |6] the calculated phonon dispersion of HgBa 2 Cu04 and Bi 2 Sr 2 Cu06 
is displayed in the main-symmetry direction A ~ (1, 0, 0), £ ~ (1, 1, 0) and A ~ (0, 0, 1). 
In addition, an interpretation of the IXS-phonon spectra for HgBa 2 Cu04 [T2l \5$\ and 
for Bi 2 Sr 2 Cu0 6 [16] is added to the figures in form of open squares and circles. 

In case of HgBa 2 Cu04 a minor deficit is visible in form of one slightly unstable 
acoustic A3 mode which is related to lattice vibrations propagating along the c-axis 
where the O z and Ba ions are vibrating in phase parallel to the CuO plane. This 
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Figure 6. Calculated phonon dispersion of E^Si^CuOg for model IT-3dim in the 
main symmetry directions A ~ (1,0,0), £ ~ (1,1,0) and A ~ (0,0,1). The symbols 
□ represent experimental results from [16] . The circle • marks the breathing modes 
at X and the half-breathing modes at A/2. For the arrangement of the panels see 
figure [3 

instability is already present in the RIM and most likely is a result of the more open 
structure of HgBa 2 Cu04 and thus of the the growing importance of non-central forces. 

The interpretation of the IXS-spectra in Refs. [12j [54] and [16], respectively, are 
based on a Voigt fit. The positions of the peaks resulting from this fit are displayed 
in the IXS-spectra shown in figure [7](a) for HgBa 2 Cu04 and [7(b) for Bi 2 Sr 2 Cu0 6 as 
open triangles. In addition, we have plotted in this figure the calculated frequencies of 
the Ax-modes in the energy range considered, including the half-breathing mode, at the 
corresponding wavevectors (full triangles). 

At first one notices that for HgBa 2 Cu04 only one phonon peak has been fitted in 
the energy range shown in figure [7(a). On the other hand, we find in this range of energy 
a phonon spectrum which is much more complex comprising up to four Ai modes, see 
figure O In the fit for Bi 2 Sr 2 Cu06 [16] only two phonon peaks have been fitted in the 
range of energy shown in figure 0(b), while according to our calculations in figure [6] there 
are up to five A± modes in this range which can interact with each other and constitute 
a highly nontrivial anticrossing scenarium. Quite general we can say that the IXS- 
spectra reveal for small ( values (q = (£,0,0)) in HgBa 2 Cu0 4 one and in Bi 2 Sr 2 Cu0 6 
two marked peaks. However, the strength of the peaks strongly decreases for larger ( 
values such that at ( ~ 0.5 the peaks are very weak and an unique interpretation seems 
not possible. Moreover, the dips in the phonon dispersion curves as extracted from the 
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Figure 7. Reconstructed experimental x-ray spectra N(E) of (a) HgBa2Cu04 for 
q = f(C,0,0) [12] and (b) Bi 2 Sr 2 Cu0 6 for q = ^((,0,0) [IB] in arbitrary units, 
respectively. The experimental spectra for the different £ are represented by the open 
and closed circles (O ). The according interpretations of the phonon frequencies 
fitted by a Voigt fit [TU [IB] are represented by open triangles (a). Our calculated 
phonon frequencies are represented by closed triangles (A). 



interpretation of the IXS-spectra at ( = 0.36 for HgBa2Cu04 [12] and at £ = 0.25 for 
Bi 2 Sr 2 Cu06 [IE] are not found in our calculations. However, as a result of the complex 
anticrossing scenario we obtain a pronounced dip at ( = 0.5 in Bi 2 Sr 2 Cu0 6 and two 
half-breathing modes below the dip with somewhat smaller frequencies, see figure [6] 
The higher one has a displacement pattern as shown in figure [3] (d) and in the lower 
one the apex oxygen and the oxygens in the ionic layer vibrate with opposite phases 
compared to figure |3] (d). This holds also true for the X-point, where two breathing 
modes occur, see figure [6] The higher one has the displacement pattern as shown in 
figure [3] (d), and in the lower one the apex oxygen and the Bismuth ion vibrate with 
opposite phases compared to figure [3] (d). 

In case of HgBa 2 Cu04 we do not find a pronounced dip for the Ai modes in the 
relevant range of energy. We obtain as a result of four interacting Ai modes two branches 
with the highest frequencies showing a downward dispersion from F to X while the 
anomalous half-breathing mode is found at a lower energy and turns out to be the 
end point of the Ai branch with the third highest frequency. Similar as in the case of 
Bi 2 Sr 2 Cu0 6 we find in HgBa 2 Cu0 4 two breathing modes at the M point. The lower 
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Figure 8. Comparison of model iT-2dim ( ) and model iT-3dim ( ). The 

three leftmost panels belong to HgBa2Cu04 and the five rightmost panels belong 
to Z^Sr^CuOg- Only irreducible representations with a different phonon dispersion 
in model II-2dim and model IT-3dim are shown. The arrangement of the panels 
from left to right according to the different irreducible representations is as follows: 
|Ai|Si|Ai||Ai|A 3 |Ei|E 3 |Ai|. 

one has a displacement pattern where the apex oxygen vibrates in opposite direction as 
compared to the higher frequency breathing mode displayed in figure figure [3] (c). 

So, we guess from our calculations for HgBa 2 Cu0 4 and Bi 2 Sr 2 Cu0 6 that the dip 
in the interpretation of the IXS-spectra may be an artefact of the complex anticrossing 
behaviour of several A x modes in the frequency range under consideration. However, we 
admit concerning our results, that our microscopic model for the electronic polarizability 
applied to the calculation of the phonon dynamics of HgBa 2 Cu04 and Bi 2 Sr 2 Cu06 is 
not fully ab initio. 

The reliability of our calculated phonon dispersion may be judged from the 
conformance of our results with the phonon branches in HgBa 2 Cu04 at lower frequencies 
measured so far and displayed in figure [5] for the model Il-3dim. 

The experimental and calculated results for the Raman and infrared active modes 
[56| Wf\ |58| |59| 160] are listed in table @] for both compounds and a good agreement is 
found for the measured and the calculated data in model Il-3dim. This fact together 
with the good results in model Il-3dim for the dispersion underlines the importance of 
the electronic polarization processes of certain ions in the ionic layers and in general of 
the physics of the third dimension perpendicular to the CuO plane. 

The effect introduced by the CF's on the O z , Hg and Bi ions outside of the CuO 
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Figure 9. Calculated phonon density of states (DOS) for (a) Hg2BaCu04 and (b) 

Bi2Sr2Cu06 for model II-3dini ( ). For comparison, the DOS of the RIM is shown 

as dotted line ( ). The symbols x and • in (a) show the experimental data [5T] at 

30 and 300K, respectively, the symbols • in (b) represent the experimental data [53] 
at 300K. All experimental data are in arbitrary units and shifted for better viewing. 



plane on the phonon dispersion by itself is demonstrated in figure [8] for both compounds 
by comparing the results of model Il-2dim (broken lines) with those of model Il-3dim 
(full lines). Only the branches are shown where large differences appear. The most 
significant renormalization of the dispersion emerges for the modes with the highest 
frequencies. 

Finally, in figure [9j the calculated phonon density of states (DOS) is shown for 
HgBa 2 Cu0 4 and Bi 2 Sr 2 Cu0 6 , respectively, and is compared with the experimental 
findings for both compounds. The theoretical results have been obtained with model 
Il-3dim. They compare quite well with the experimental results for both, the width of 
the spectra and their peak structure. Of course the assignment of the peaks is difficult 
in detail because the experimental data are broadened. 

A calculation of the partial phonon density of states, not shown in this work, 
demonstrates that the DOS in the lower frequency range is dominated by Cu-, Ba-, 
Hg-vibrations and by Cu-, Bi-, Sr-vibrations, respectively. On the other hand, the high- 
frequency part of the spectrum is nearly exclusively governed by the oxygen modes. 

From figure [9] (a) and (b) also the difference of the DOS in the model Il-3dim and 
the RIM generated by the electronic polarization processes can be extracted. From the 
calculated results we find a significant decrease of the width of the spectra together with 
a characteristic redistribution of spectral weight particularly in the high-frequency part. 
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The latter is dominated by strongly coupling high-frequency oxgen vibrations living in 
this part of the spectra and suffering a corresponding softening compared to the RIM. 
The pronounced peaks which are generated by this softening also mean that there is 
considerable phase space for this strongly coupling modes. 

4. Summary and conclusions 

We have shown that for a reliable description of phonon dynamics and charge response 
in HgBa2Cu04 and Bi 2 Sr 2 Cu06 a broad set of orbital degrees of freedom (Cu3d, 
Cu4s, 02p) must be taken into account to model the electronic structure of the CuO 
plane. Moreover, for an ample representation of the real three-dimensional compounds 
electronic degrees of freedom beyond the strictly two-dimensional electronic structure 
must be considered. 

We found that in the link between the CuO plane and the ionic layers of the 
cuprates electronic polarization processes of CF type are essential at certain ions. In 
HgBa 2 Cu0 4 O z 2p and Hg6s CF's and in Bi 2 Sr 2 Cu0 6 O z 2p and Bi6p CF's, respectively, 
are significant processes. CF's on these ions are promoted by the short distance of the 
2 -Hg and O^-Bi bond. Allowance of CF's on these ions outside of the CuO plane 
(model Il-3dim) leads to a much better depiction of the phonon dispersion curves, the 
Raman active and the infrared active modes, in particular for those where the apex 
oxygen vibrates along its bond with the Hg or the Bi ion, respectively. 

As far as the electronic degrees of freedom in the CuO plane are concerned we 
have pointed out that besides the localized Cu3rf state and the 02p state also the more 
delocalized Cu4s state is needed in the charge response for an adequate representation 
of the large softening of the high-frequency oxygen bond-stretching modes and the 
strong coupling to the electrons. Such an effect also has been found in our previous 
calculations for La 2 Cu0 4 and Nd 2 Cu0 4 . The strong nonlocal coupling and the related 
strong renormalization of these modes seems to be generic in the cuprates. 

Altogether, our calculated results for model Il-3dim are partly in good agreement 
with the phonon dispersion measured so far in HgBa 2 Cu04 and Bi 2 Sr 2 Cu06 and also 
with the Raman active and infrared active modes. 

However, we have pointed out that an interpretation of the IXS measurements 
of HgBa 2 Cu0 4 and Bi 2 Sr 2 Cu0 6 in the high-frequency part of the spectrum and the 
assignment of dips in the corresponding phonon dispersion curves cannot be supported 
by our calculations. This is, because our investigations exhibit that in the frequency 
range under study we have multiple interacting phonon modes of A x symmetry (four in 
HgBa 2 Cu0 4 and five in Bi 2 Sr 2 CuOe) which create a complex anticrossing scenario that 
possibly cannot be studied correctly within a fit using just one or two modes. 

Finally, we have calculated the DOS of both compounds in satisfactory agreement 
with experiment. We found a significant redistribution of spectral weight due to the 
electronic polarization processes in particular in the high-frequency part of the spectrum. 
This energy range is dominated by the strongly coupling oxygen modes signaling a 
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considerable amount of phase space for strong electron phonon coupling. 
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